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Abstract 

Small superconducting grains are discussed in the frameworks of both the reduced 
BCS Hamiltonian and the Universal Hamiltonian. It is shown that fluctuations of 
electrons in levels far from the Fermi energy dominate superconducting properties 
in small and ultrasmall grains. Experimental consequences related to the spin sus- 
ceptibility and persistent currents of grains and rings with weak electron-electron 
interactions are discussed. 

Key words: superconductivity, granular, fluctuations 
PACS: 74.81.Bd, 74.25.Bt, 73.23.Ra, 74.20.Fg 



1 Introduction 



The paradigmatic framework in the study of superconducting fluctuations in 
small grainspQ is the reduced BCS Hamiltonian 

3,<r=± i,3 



Here A is the dimensionless interaction constant, d is the mean level spacing, 
and the indices i,j correspond to doubly degenerate time reversed states of 
energy Ep — ujd < e < E-p + up,. 

While the mean field (BCS)[2J solution of this Hamiltonian was extremely 
powerful in explaining the properties of bulk superconductors, grains of fi- 
nite size pose additional challenges. In such grains the number of electrons is 
fixed, their size may be smaller than the coherence length £ and the magnetic 
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penetration depth, and zero resistance is not achievable. Still, already in 1959 
Anderson postulated that superconducting grains maintain superconducting 
properties for sizes much smaller than £, down to the size in which the single 
electron level spacing d equals the bulk energy gap A (the " Anderson size" [3] ) • 
For clean grains this size is of order £ 1//3 Aj!/ 3 <C £. 

Much of the recent theoretical interest in the study of small superconducting 
grains (see Ref.pQ and references therein) was initiated by the experimental 
work of Ralph Black and Tinkham[l], who found a parity dependent gap in 
the excitation spectrum of single "small" Al grains (with d < A), but not for 
"ultrasmall" grains having d > A. More recently it was also shown that the 
Meissner effect sharply disappears at grain sizes consistent with the Anderson 
size [51161171 . 

While superconducting characteristics reminiscent of bulk properties indeed 
vanish at the Anderson size, signatures of pairing correlations persist to smaller 
sizes. These signatures are a result of superconducting fluctuations of electrons 
at energies larger than the gap energy, up to the cutoff energy given by c^d- 
For example, the spin susceptibility of ultrasmall superconducting grains is 
predicted to exhibit a re-entrant behavior both as function of temperature [8] 
and as function of magnetic field j9|T0| . with a long tail, persisting up to the 
temperature/magnetic field equivalent of u D . The superconducting fluctua- 
tions of electrons further than A from the Fermi energy Ep affect significantly 
also grains in an intermediate regime, where A 2 /cg»d < d < A. In particu- 
lar, in this regime the energy gain of the system by the attractive interaction 
is much larger than that given by the mean field treatment of BCS [9], and 
the same is true for the pairing parameter [llj. Moreover, it was shown that 
these far level fluctuations affect differently single particle and collective prop- 
erties, and therefore one has to define correspondingly two different pairing 
parameters pTj when discussing small superconducting grains. 

While the above properties, and specifically their dependence on the high en- 
ergy cutoff, stem from the exact solution of the reduced BCS hamirtonian |12|[T3] 
the question to their relevance to real superconducting grains may arise [14j. 
This question is of particular relevance since specific predictions for the ex- 
perimental consequences of the fluctuations of the far levels (e.g. the long tail 
of the spin susceptibility discussed above) are made. Here we discuss the sig- 
nificance of the superconducting fluctuations of the far levels by considering 
both the reduced BCS model and the model of the Universal hamiltonian[15j. 
While the reduced BCS Hamiltonian is an effective Hamiltonian whose va- 
lidity in describing high energy properties may be questioned, the Universal 
Hamiltonian was shown, using renormalization group approach, to control the 
low energy physics of metallic grains with weak interactions and large dimen- 
sionless conductance [T6l[TT] . We thus establish the significance of the super- 
conducting fluctuations of the far levels up to the high energy cutoff of the 
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Universal Hamiltonian given by the Thouless energy Efh- We then discuss the 
significance of the Debye energy in general, and for the problem of persistent 
currents in particular. 



2 Reduced BCS Hamiltonian and Universal Hamiltonian - exact 
solution 

The reduced BCS Hamiltonian (JTJ), with finite number of electrons, was solved 
exactly by Richardson [T2iri3] . The structure of the solution takes advantage 
of the fact that the interaction scatters only pairs and not single electrons. 
Thus, the solution is obtained separately for subspaces of the Hilbert space 
defined by the identity of the singly occupied levels. The singly occupied levels 
are neglected at first, the problem of N pairs in M states is solved, and the 
singly occupied levels are then trivially added to the solution (see Refs. [T2l[T3"] 
as well as Ref.[18] for details). 

While the reduced BCS Hamiltonian is motivated by the specific phonon- 
mediated attractive electron-electron interaction, it was recently shown that 
any metallic grains with large dimensionless conductance g = E^h/d, weak 
interactions, and negligible spin-orbit interaction, can be described by the 
Universal Hamiltonian [15], which includes only three interaction parameters 

H = E en4,«rCn,* + E C N 2 + J C T*T + J S S 2 . (2) 

n,cr 

Here N = £ no . 4 j0 .c^ )(7 is the number operator, S = | En v 4, ff V c «y is 
the total spin operator, and T = J2n c n ,-C n ,+ is the pair annihilation operator. 
The index n spans a shell of doubly degenerate time reversed states of energy 
E-p — E-rh < £n < Ep + E^. E c is the charging energy and J c ( s ) = X c ( s )d, where 
A c and X s are the dimensionless interaction parameters in the Cooper channel 
and in the spin channel respectively. 

Interestingly, for isolated grains the solution of the Universal Hamiltonian is 
given by Richardson's solution for the reduced BCS Hamiltonian. First, the 
Coulomb term can be neglected, since the number of electrons in the iso- 
lated grain remains constant. Then, other than a different cutoff energy (Ej^ 
compared to Wj)), one remains with the reduced BCS Hamiltonian with the 
additional exchange term. The latter commutes with the rest of the Hamilto- 
nian. Moreover, the pairing interaction involves solely the paired levels, while 
the spin interaction involves solely the singly occupied levels. One can then 
obtain the solution for the Universal Hamiltonian by following the steps in 
Richardson's solution of the reduced BCS Hamiltonian, with the additional 
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consideration of the spin term for the singly occupied levels. 



3 Large contribution of the far levels 

Considering the reduced BCS Hamiltonian (p^), and using Richardson's exact 
solution, the condensation energy of small metallic grains was calculated as 
function of the coupling constant A [9j. The condensation energy was defined 
as the difference between the energy of the Fermi state and the real ground 
state of the system, i.e. 

-Econd(A) = -E'F.g.s(A) — E gs _(X). (3) 

For bulk superconductors this energy is given by A 2 / (2d), is extensive, and 
is a non-analytic function of A. For finite size grains it was found that the 
condensation energy is analytic at A = 0, and is very well estimated by -E C ond — 
A 2 / (2d) +ln 2A 2 co>d- The first term gives the contribution of the levels within A 
of Ep, and the second, perturbative term, is the contribution of the far levels. 
Interestingly, due to unique analytical properties of the condensation energy, 
the perturbative term expresses the contribution of the far levels not only 
within the regime of validity of perturbation theory (A < 1/ In N, or d > A) , 
but also in the regime where d < A [9]. An immediate outcome of this result is 
that for a large intermediate regime, A 2 /uo^ < d < A, the condensation energy 
is much larger than that given by the BCS term. Note that the perturbative 
term is intensive, and therefore negligible as the size of the grain becomes 
large. 

Similarly, the contribution of the far levels is found to be significant when 
considering generalizations of the bulk order parameter that are suitable for 
finite size grains [TT] (see also Ref. pQ and references therein). For all standard 
definitions of the order parameter the far level contribution results in a term 
linear in u-q, and the order parameter being much larger than its mean field 
value in the intermediate regime defined above (see Ref. [TT] for details). Thus, 
in contrast to the common belief that the order parameter turns from being 
extensive to being intensive at the Anderson size (d ~ A), it actually becomes 
intensive already at a much larger size, i.e. when d A 2 /ujd. Furthermore, 
the significant contribution of the far levels separates collective properties 
of the superconductor such as the condensation energy, from single particle 
properties such as the energy gap for single particle excitations. While in bulk 
superconductors both energies are related to the bulk order parameter A, 
for small grains one has to define two different order parameters to describe 
collective and single particle properties[UJ. While the cutoff energy affects the 
former linearly, it affects single particle properties only logarithmically. 
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Interestingly, using the similarity of the reduced BCS and the Universal Hamil- 
tonians, and their exact solution, the results discussed above for the reduced 
BCS Hamiltonain are immediately applicable to the Universal Hamiltonian, 
with the only change of the cutoff from to E^ h . Thus, as long as the grains 
obey the conditions of applicability of the Universal Hamiltonian described 
above, one obtains the intermediate regime, defined now by A 2 /i?Th < d < A, 
in which e.g. the condensation energy is much larger than the mean field BCS 
value, and therefore intrinsic. The predictions of re-entrant spin susceptibility 
as function of temperature [S] and magnetic Held !), 10 . which result from the 
pairing correlations of the far level, can also be made on the basis of the va- 
lidity of the Universal Hamiltonian. Importantly, all the above characteristics 
do not depend on the interaction being constant up to the upper cutoff, nor 
do they depend on the existence of a sharp cutoff. It is sufficient that one 
can bound the interaction from below by CiA c within a window of c^E^ from 
E F , where Ci,c 2 are constants of order unity. Thus, based on the validity of 
the Universal Hamiltonian one can argue for the applicability of the above 
results for the condensation energy and spin susceptibility for real physical 
systems. Note, that within the regime of applicability of the Universal Hamil- 
tonian one can uniquely relate the excess spin susceptibility as function of the 
magnetic field to the existence of pairing correlations, as the exchange term 
can not account for such behavior. Furthermore, by examining the behavior 
of the spin susceptibility of an ensemble of small metallic grains the interac- 
tion parameters of the Universal Hamiltonian can be determined |10j. Another 
consequence of the above analysis is that the interesting question of whether 
the noble metals have a weak attractive interaction and therefore are super- 
conductors albeit with a very low T c can be addressed as well, in small grains 
and by extrapolation in bulk. This is since unlike T c which is exponentially 
small in A c , the excess spin susceptibility as a function of the magnetic field 
is quadratic in A c and therefore detectable for small A c . 



4 Persistent currents and the significance of c^d 

In the discussion above it is argued that the interesting physics related to the 
pairing correlations of the far levels is independent from the validity of the 
reduced BCS model in describing the properties of small grains at energies of 
the order of uj d . Still, the energy scale of u D is a physical energy scale, related 
to the retardation of the phonon mediated interaction. Especially since, as was 
mentioned above for the Universal Hamiltonian, a constant interaction and a 
sharp cutoff at cud are not necessary for the applicability of the above results, 
it is plausible that indeed the energy scale of wd dictates the magnitude of the 
condensation energy of small grains, the boundary of the intermediate regime, 
and the extent of the tail in the excess spin susceptibility. 
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However, in relation to persistent currents, a rigorous understanding of the in- 
teraction is even more crucial. Since the persistent current in a normal metal 
ring is given by the derivative of the energy with respect to the flux, it is plau- 
sible that for metals with weak attractive interaction, the large contribution of 
the far levels to the condensation energy will affect the persistent current, and 
will result, within the reduced BCS model, in a large persistent current related 
to the energy scale of uj^. Indeed, in Ref. [19] it was shown that the value of 
the derivative of the persistent current at zero flux is much larger within the 
reduced BCS model in comparison to the value obtained within the standard 
theory of momentum independent interaction [20] . Here, however, the differ- 
ence in magnitude of the persistent current is related to the high energy cutoff 
being un rather than E^, and thus is crucially related to the specifics of the 
model considered, and to the exact form of the interaction. In particular, the 
BCS interaction assumes that only time reversed states interact, i.e. that the 
total incoming momentum of the scattered electrons q must be zero. While 
this is a simplified form of the momentum dependence of the interaction, the 
subsequent result of the much larger magnetic reponse[19] in comparison to 
the value obtained within the momentum independent picture points to the 
importance of taking correctly the dependence of the interaction on q, espe- 
cially since a significant q dependence of the attractive interaction is motivated 
by the retardation of the phonon mediated interaction [191121] . We believe that 
a rigorous understanding of the form of the interaction, and in particular its 
q dependence, could lead to a better understanding of the long standing ques- 
tion regarding the value of the ensemble averaged persistent current in small 
metallic rings. 

This work was supported by NSERC of Canada and PITP. 
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